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1.  INTRODUCTION 


This  is  SARA’s  3rd  Quarterly  Report  for  “Breakthroughs  in  Low-profile  Leaky-Wave  HPM  Antennas,”  a 
37-month  Basic  Research  effort,  sponsored  by  the  US  Office  of  Naval  Research  (ONR).  This  work 
includes  fundamental  theoretical  analyses,  numerical  modeling,  and  related  basic  research.  Objectives 
include  to  discover,  identify,  investigate,  characterize,  quantify,  and  document  the  performance,  behavior, 
and  design  of  innovative  High  Power  Microwave  (HPM,  GW-class)  antennas  of  the  forward-traveling, 
fast-wave,  leaky-wave  class. 


1 .1 .  Overview  of  Previous  Activities  (1st  and  2nd  Quarter) 

During  the  first  quarter,  we  prepared  and  established  useful  equations  and  algorithms  for  predicting 
reflections  and  transmission  of  incident  TE  waves  from  parallel-wire  grills,  dielectric  windows,  and 
combinations  of  wire  grills  with  dielectric  windows,  in  problems  reducible  to  purely  H-plane  (2D) 
representations.  We  then  applied  this  theory  to  guide  the  design  of  high-gain  configurations  (again, 
limited  to  2D,  H-plane  representations)  for  linear,  forward  traveling-wave,  leaky- wave  antennas.  The 
theory  built  upon  equivalent  circuit  methods  and  wave  matrix  theory,  which  provided  useful  formalisms 
upon  which  we  continue  to  build. 

During  the  second  quarter,  we  pursued  initial  extensions  of  the  previous  work  into  three  dimensions,  in 
order  to  include  phenomena  with  E-plane  dependencies.  We  succeeded  in  adding  into  the  wave-matrix 
formalism  the  reflection/transmission  properties  associated  with  the  transition  to  free  space  from  &  finite- 
width  leaky-wave  channel,  including  the  edge-tapering  essential  to  HPM  applications.  These  geometric 
aspects  do  not  arise  in  analyses  confined  to  the  H-plane  alone.  As  expected,  our  3D  analyses  were 
somewhat  more  reliant  on  numerical  models  than  in  the  2D  analyses,  due  to  the  greater  complexity  of 
identifying  and/or  building  practical  analytic  approaches  capable  of  addressing  true  3D  geometries  of 
interest.  For  more  information,  we  refer  the  reader  to  our  Quarterly  Report  # 1  and  Quarterly  Report  #2. 

1 .2.  Overview  of  Recent  Activities  (3rd  Quarter) 

During  the  third  quarter,  we  pursued  a  variety  of  continuing  and  new  analyses/design  activities,  many  of 
which  are  still  in-progress.  One  emphasis  during  this  period  was  exploration  of  channel-to-channel 
coupling  (aka,  mutual  coupling)  which  (as  we  have  noted  earlier)  is  an  important  design  concern,  since  it 
can  impact  antenna  performance  significantly  in  terms  of  gain,  peak  power-handling,  and  impedance 
matching.  Our  approach  here  leverages  mostly  numerical  methods,  but  with  intuitive  considerations. 

This  report  also  includes  a  pleasingly-simple  derivation  (in  2D  geometry)  of  an  equation  that  tells  how  to 
position  the  leaky-grill  relative  to  the  opposite  wall  to  compensate  for  changes  in  the  phase  velocity  due 
to  leakage  through  the  grill.  In  the  form  derived,  the  result  appears  to  be  generalizable  to  3D  geometries. 

In  addition  to  those  topics,  we  discuss  some  of  the  aperture  geometries  (shapes  and  curvatures)  currently 
being  investigated.  Section  3  describes  the  technical  work  performed  this  quarter  in  more  detail. 

2.  STATUS  OF  THE  PLAN/SCHEDULE  AND  FUNDING 

Figure  1  below  maps  out  the  overall  program  plan,  for  quick  reference.  As  of  the  time  of  this  report,  we 
are  continuing  Tasks  2.1,  2.2,  3.1,  and  are  initiating  Tasks  2.3  and  3.2.  The  latter  have  been  somewhat 
delayed  relative  to  our  initial  plans,  but  we  plan/expect  to  catch  up  over  the  next  quarter. 

The  subject  contract  was  awarded  on  9/18/2013  and  has  an  end  date  of  10/17/2016.  The  total  contract 
value  is  $868,350,  with  current  (per  P00003  signed  on  4/24/2014)  allotted  funding  of  $406,530. 
According  to  SARA’s  accounting  system  as  of  June  13,  2014,  expenses  (including  fee)  have  totaled 
$211,116,  thus  leaving  $195,414  remaining. 
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Plan  of  Action  and  Milestones  (POA&M)  (updated  June  2014) 


Activity  Name 
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♦  ♦ 

10/14, 10/15 

J  ♦ 

♦ 

-10/4/16 

♦ 

-10/1 

Ll/16 

♦ 

1.0  Program  Management  &  Reporting 

1.1  Kickoff  Meeting  (1) 

1.2  Quarterly  Reports  (11  tot) 

1.3  Annual  Review  Meetings  (2  tot) 

1.4  Final  Report  (1) 

1.5  Final  Review  Meeting  (1) 


2.0  Prep  Fundamental  Analyses  &  Models 

09/18/13 

02/28/16 

2.1  Establish/document  FAWSEA  theory 

09/18/13 

06/30/14 

2.2  Add/generalize  theory  to  include  CAWSEA 

11/01/13 

09/30/14 

2.3  Extend  theory  to  include  AAWSEA 

04/01/14 

03/31/15 

2.4  Add  RAWSEA-specific  considerations 

10/01/14 

06/30/15 

2.5  Extend  theory  to  new  designs  (see  below) 

04/01/15 

02/28/16 

3.0  Establish  Optimal/Recommended  Designs 

01/01/14 

06/30/16 

3.1  FAWSEA-based 

01/01/14 

09/30/14 

3.2  CAWSEA-based 

04/01/14 

12/30/14 

3.3  AAWSEA-based 

10/01/14 

06/30/15 

3.4  RAWSEA-based 

04/01/15 

12/31/15 

3.5  Designs  leveraging  new  features  (see  below) 

10/01/15 

06/30/16 

4.0  Develop  &  Document  New  Designs 

07/01/14 

08/01/16 

4.1  New  variant  #1  (e.g.,  "Pinched..."  ->  PAWSEA) 

07/01/14 

03/30/15 

4.2  New  variant  M2 

04/01/15 

12/31/15 

4.3  New  variant  M3  (+ any  others) 

01/01/16 

08/01/16 

Figure  I,  Overall  program  plan 


3.  RESEARCH  PERFORMED  THIS  PERIOD 


3.1 .  Investigations  of  Mutual  Coupling 

It  would,  of  course,  be  very  helpful  if  the  design  of  forward-traveling  leaky-wave  HPM-capable  antennas 
could  proceed  without  consideration  of  the  effects  of  mutual  coupling  between  channels.  However,  such 
effects  can  and  do  have  significant  impact  on  resulting  performance,  leading  to  increased  design  time, 
uncertainty,  and  constraints  on  design  choices.  Mutual  coupling  has  long  received  attention1  in  the 
literature,  but  the  subject  remains  theoretically  challenging,  although  simple  parallel  dipole  elements  have 
been  analyzed  theoretically  in  much  detail1.  The  papers  by  S.  Nishida2  are  perhaps  the  most  relevant  to 
our  own  work,  but  do  not  include  essential  geometric  features  (such  as  rounded  edges,  and  not  to  even 
mention  dielectric  windows)  required  for  HPM-capable  designs. 

In  analogy  to  the  numerical  study  investigating  the  impact  of  rounded  edges  in  the  E-plane  as  described  in 
our  previous  report,  we  performed  another  2D  FE-based  numerical  study  investigating  mutual  coupling 
between  radiating  openings  with  rounded  (aka,  filleted)  edges.  This  study  investigated  the  dependence  of 
mutual  coupling  upon  channel  separations,  choice  of  fillet  radius  at  the  channel  openings,  and  frequency. 
Figure  2  and  Figure  3  summarize  the  setup  of  the  2D  model.  Intermediate  boundaries  to  gently  scale  the 
mesh  and  a  perfectly  matched  layer  (PMF)  outer  boundary  were  incorporated  to  improve  speed  and 
accuracy,  respectively.  One  hundred  geometries  were  considered  at  33  frequency  values,  for  a  total  of 
3,300  runs.  A  plot  showing  all  3,330  “raw”  values  of  the  computed  magnitude  of  S2i  is  provided  in 
Figure  4,  while  projections  of  these  values  against  some  parameters  of  interest  are  shown  in  Figure  5. 

The  closest  to  a  “single-parameter”  dependence  of  S2i(dB)  is  seen  for  the  channel  separation/wavelength. 


1  E.g.,  see  Oliner,  A.A.  and  R.G.  Malech,  “Radiating  Elements  and  Mutual  Coupling,”  “Mutual  Coupling  in  Infinite 
Scanning  Arrays,”  and  “Mutual  Coupling  in  Finite  Scanning  Arrays,”  —  Chaps.  2,  3,  and  4  respectively  of  Array 
Theory  and  Practice ,  Vol.  II  of  Microwave  Scanning  Antennas ,  Ed.  by  R.C.  Hansen,  Peninsula  Publishing,  Los 
Altos,  CA,  1985. 

2  Nishida,  S.,  “Coupled  Leaky  Waveguides  I:  Two  Parallel  Slits  in  a  Plane,”  IRE  Trans.  Ant.  and  Propagat .,  May, 
1960,  pp.  323-330.  See  also  Nishida,  S.,  “Coupled  Leaky  Waveguides  II:  Two  Parallel  Slits  in  a  Cylinder,”  IRE 
Trans.  Ant.  and  Propagat .,  July,  1960,  pp.  354-360. 
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Parameter 
Chan,  width: 

Rfil: 

Gap: 

Frequency: 


Range 

Fixed  at  5.0  cm 

From  1.0  to  10.0  cm,  in  1.0  cm  steps 
From  1.0  to  10.0  cm,  in  1.0  cm  steps 
50  to  200  MHz  in  25  MHz  steps,  plus 
250  to  1500  MHz  in  50  MHz  steps 


#  of  cases 

1 

10  \  . 

10  >  -7  100  geometries 

33 

3,300  2D-model  runs 


2D  Parameterization: 


PEC  symmetry  / me 

-0.25  -0.2  -0.15  -0.1  -0.05  0  0.05  0.1 


Chan,  sep  =  Chan,  width  +  2Rfjl  +  Gap 


Extremes  of  the  geometry: 


Rfii  =  10cm 


(Mi  Rfi|  - 1  cm 
°”  Gap  =  10cm 

4*1 

•J 


01 


0 

-44  53  -41 

Rfii  = 10cm 

0  41  42 

Gap  =  10cm 

1 

ft 

it 

A 

44  4)  4]  41  0  01  02 


Figure  2.  Parameter  definitions  in  2D  E-plane  TEM-wave  mutual-coupling  study. 


Mesh  (example) 


Sample  run 

freq(l)=10e8  Surface:  Magnetic  field,  z  component  (A/m) 


-0.5 


2.2 

2 

1.8' 

1.6 

1.4 

1.2' 

r 

0.8' 

0.6 

0.4' 

0.2' 

O' 


2.2 

2 

1.8 

1.6 

1.4 

1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 


▲  0.501 
0.5 
0.4 
0.3 
0.2 
0.1 
0 

-0.1 

-0.2 

-0.3 

-0.4 

-0.5 
▼  -0.501 


From:  Guides_coupling_2D_testbed3_onecase  .mph 

Figure  3.  Computation  space  employed  in  2D  E-plane  TEM-wave  mutual-coupling  study. 
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S21,<dB>  ,  .  ,S21,<dB>,  ,  ,  ,  Increasing  mutual  coupling 


Most 

Global:  Adjstd  S21dB  (dB) 


Rfii  ~  1.0cm, 
Gap  =  1.0cm 


-12 


-13 


-14 


-15 


-16 


-17 


-18 


-19 


-20 


-21 


-22 


-23 


-24 


-25 


freq  From:  Guidcs_coupling_2D_lcslbed3.mph 

Among  the  above  cases,  the  strongest  mutual  coupling  occurred  at 
Rfii  =  1.0cm,  Gap  =  1.0cm,  and  f  =  850  MHz  ->  S21  =  -11.436  dB 


•  As  expected,  strongest 
mutual  coupling  occurs 
when  the  channels  are 
closest  together;  weakest 
coupling  when  they  are 
farthest  apart. 

•  Coupling  becomes  weak 
at  very  low  frequencies 
due  to  the  unavoidably 
poor  impedance  match  of 
each  channel  to  its  own 
aperture. 

•  Coupling  also  weakens 
at  sufficiently  high 
frequencies  due  to 
greater  effective  channel 
separation,  rel.toA,. 

•  Ripples  in  coupling 
curves  vs.  freq  are  seen 
in  the  vicinity  of  channel 
separations  ~  nX/2. 


Figure  4.  Computed  S2i  (dB)  vs,  frequency  (all  3,300  cases) 
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Gap/wavelength 
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Figure  5.  Projections  of  values  in  Figure  4  vs.  various  parameters  scaled  by  wavelength. 
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Application  to  3D 

To  better  understand  the  implications  of  the  2D  model  in  a  3D  context,  we  pursued  the  following 
approach  with  the  help  of  3D  models: 

•  Based  on  the  2D  analyses,  consider  a  3D  FAWSEA  design  likely  to  exhibit  strong  mutual  coupling, 
so  that  any  problematic  effects  of  mutual  coupling  will  be  easy  to  observe. 

•  In  particular,  noting  the  highest  value  of  S2i(dB)  in  Figure  4,  consider  a  FAWSEA  with  fTEM  =  850 
MHz,  with  Rffl  =  1.0cm,  Gap  =  1.0cm,  and  channel  width  =  5cm  (yields  channel  spacing  =  8cm). 

•  Choose  the  frequency  in  the  3D  model  to  account  for  the  tilted-wave  incidence  (vs.  normal  incidence 
in  the  2D  TEM  model).  I.e.,  ^adj  =  A3cos((|)in),  where  A,’  comes  from  the  2D-model’s  fTEM=850  MHz. 
If  we  choose  (J)in  =  30°  (a  typical  value,  as  noted  previously),  we  obtain  f3D  =  981.5  MHz. 

•  Apply  the  3D  FAWSEA  design  algorithms  &  scripts  from  our  previous  work,  but  without  a  dielectric 
window  (for  simplicity),  to  set  the  diameters  of  the  grill  wires. 

•  Compare  single-channel  vs.  three-channel  3D  configurations,  the  latter  variously  with: 

(a)  all  three  channels  driven, 

(b)  only  a  single  side-channel  driven,  and 

(c)  only  the  center  channel  driven. 

We  selected  the  following  as  3D  FAWSEA  design  parameters  (i.e.,  the  inputs  to  our  FAWSEA  design 
tools/scripts3)  :  foperating  =  981.5  MHz,  4>in  =  30°,  chan,  length  =  2.0  m,  chan,  width  =  5.0  cm,  Rm  =  1  cm, 
Gap  =  1  cm  (->  chan,  separation  =  8  cm),  and  a  constant  grill- wire  spacing  =  6cm.  Figure  6  shows  the 
computed  wire  diameters: 


*Applied  routines:  fawsea_script_3D0_radiused.m\N\th  myRpowF.m 


Position 

Diameter 

Position 

Diameter  (cont 

0.06 

0.00708204 

0.96 

0.00513925 

0.12 

0.00701018 

1.02 

0.00496521 

0.18 

0.00690004 

1.08 

0.00475043 

0.24 

0.00679323 

1.14 

0.00455415 

0.3 

0.0066613 

1.2 

0.00431537 

0.36 

0.00654014 

1.26 

0.00410282 

0.42 

0.00643618 

1.32 

0.00387251 

0.48 

0.00630619 

1.38 

0.00361951 

0.54 

0.00620094 

1.44 

0.00334765 

0.6 

0.00604983 

1.5 

0.00303933 

0.66 

0.00591977 

1.56 

0.0027051 

0.72 

0.00577664 

1.62 

0.00237408 

0.78 

0.00562465 

1.68 

0.00198597 

0.84 

0.00547871 

1.74 

0.00156414 

0.9 

0.00530538 

1.8 

0.00108827 

1.86 

0.000577093 

With  channel  depth  (initial)  =  0.17634777  m. 

(Reduced  by  2.25%  later,  to  yield  clean  30°  beam  tilt.) 


Figure  6.  Computed  wire  diameters  vs  position  for  the  3D  FAWSEA  (with  no  window) 
selected  for  initial  use  in  our  mutual-coupling  study. 


Figure  7  shows  a  3D  model  of  a  single-channel  version  of  this  antenna,  including  an  absorbing 
termination4.  Adding  two  closely-adjacent  channels  yields  the  model/results  shown  in  Figure  8.  If  there 
was  no  interaction  between  the  channels,  we  would  expect  the  gain  to  increase  by  10*logi0(3)  =  4.77dB. 
But  instead,  we  find  the  gain  has  increased  by  only  3.23dB,  the  H-plane  peak  direction  has  shifted  to 
28.8°  instead  of  30°,  and  a  slightly  smaller  fraction  of  the  source  power  is  actually  being  radiated.  Thus, 
mutual  coupling  in  this  example  has  had  a  deleterious  effect. 


3  Applied  scripts/functions:  fawsea_script_3D0_radiused.m  with  myRpowF.m 

4  In  practical  designs,  the  termination  would  be  a  conductor,  but  an  absorbing  termination  allows  us  to  better 
separate  mutual  coupling  effects  from  effects  due  to  reflections  of  non-radiated  waves  reaching  the  terminations. 
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freq(l)=9.815e8  Slice:  Electric  field,  x  component  (V/m)  Slice:  Electric  field, 

Slice  vipw 
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From  3D  FAWSEA_coupling_one_isolated  June5  2014  mph 
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|  Ex  |  along  center  of  the  channel 


Far  Field:  Gain  Pattern  in  H-plane  (dB) 


•80  -60  -40  -20  0  20  40  60  80 

(|vangle  (deg) 


This  is  a  reasonable  single-channel  FAVJSEk 
example  (but  with  no  window). 

Note  (for  future  reference)  that  G  =  13.244  dB 
at  30.0°  for  this  isolated,  single-channel,  case. 


Figure  7.  Case  of  one  isolated  3D  FAWSEA  channel  (opening  onto  a  conducting  plane). 
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Far  Field:  Gain  Pattern  In  H-plane  (d8) 


Observed  departures  from  ideal  arraying: 

•  Center  chan  amplitudes  decreased;  outer  chan, 
amplitudes  increased,  rel.  to  an  isolated  chan. 

•  H-plane  gain  peak  shifted  to  28.8°  (was  30.0°) 

•  Gain  increased  by  only  +3.23  dB  (vs.  +4.77  dB, 
which  would  occur  if  ideal  combining)  relative 
to  single-channel  configuration. 

•  Slightly-reduced  fraction  of  Psrc  radiated. 


Figure  8.  Three  tightly-spaced  channels,  all  driven.  (Channel  spacing  =  8  cm) 
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One  perspective  by  which  we  can  partially  understand  the  reduction  in  gain  is  as  follows:  First,  recall  that 
the  single  isolated  channel  exhibited  gain  =  13.244  dB  (GN  =  21.106).  If  we  associate  this  with  an 
effective  area,  we  obtain  Aeff  =  GN; i2 /4tt  =  15.67  cm2.  Since  the  channel  is  2m  long,  this  implies  an 

effective  width  Wejj -  7.835  cm.  However,  that  ignores  the  efficiency  factor  associated  with  the  beam  tilt 
angle.  And  if  we  include  that  term,  the  effective  length  is  reduced  to  2m*cos(30°)  and  the  effective  width 
is  increased  to  Wejy  =  9.047  cm.  This  width  is  larger  than  the  spacing  between  the  channels  in  the  model 
of  Figure  8.  So  it  makes  sense  that  this  three-channel  configuration  is  just  too  closely -packed,  i.e.,  there 
is  insufficient  aperture  area  available  (at  least,  for  the  middle  channel)  to  deliver  the  full  array  gain. 
Further  understanding  follows  from  the  models  in  Figure  9  and  Figure  10,  which  explore  the  behavior 
when  only  a  side-channel  is  driven  and  when  only  the  center-channel  is  driven,  respectively.  From  Figure 
9,  we  see  that  the  presence  of  the  two  other  channels  has  a  moderate  impact  on  the  performance  of  the 
driven  side-channel  (e.g.,  the  gain  falls  to  12.86  dB).  But  if  just  the  center  channel  is  driven  (Figure  10), 
the  channels  immediately  on  either  side  impact  it  severely,  with  the  gain  falling  to  9.04  dB. 


freq(l)=9.815e8  Slice:  Electric  field,  x  component  (V/m)  0 

Slice:  Electric  field,  x  component  (V/m) 

Slice  view,  Ex  _ , 


(driven  side-chan  mid-plane) 


Far  Field:  Gain  Pattern  In  H-plane  (dB) 


Figure  9.  Three  tightly-spaced  channels, 


Concerns  compared  to  single  isolated  channel: 

•  Coupling  is  evident  to  neighboring  channels. 

•  H-plane  gain  peak  shifted  to  29.2°  (was  30.0°) 

•  Gain  fell  to  12.86  dB  (vs.  13.244  dB  for  an 
isolated  single-channel  configuration). 

•  Only  78%  of  source  power  radiated.  [22%  of  Psrc 
was  variously  reflected  or  lost  into  the  termina¬ 
tions  of  the  driven  and  other  two  channels.] 

!  side-channel  driven.  Chan,  spacing:  8  cm. 


To  help  clarify  what  to  do,  Figure  1 1  revisits  some  of  the  curves  from  Figure  4,  highlighting  a  few  of  the 
RfU  =  1.0  cm  cases  with  different  gaps.  Increasing  the  gap  from  1cm  to  3cm  (which  means  the  channel 
separation  increases  from  8cm  to  10cm)  should  reduce  the  mutual  coupling  by  ~ldB.  Is  that  enough?  No. 
Figure  12  repeats  the  calculation  in  Figure  10,  but  for  this  modestly -increased  spacing.  A  collapse  in  gain 
for  the  center  channel  is  still  there,  albeit  less  dramatic.  Clearly,  the  channels  are  still  too  close  together. 
Figure  13  shows  the  same  calculation  for  a  channel-to-channel  array  spacing  of  12cm  and  Figure  14 
shows  it  for  a  14cm  spacing.  This  last  case  might  actually  seems  pretty  good  at  first  look,  since  the  gain 
of  the  center  channel  is  very  nearly  that  of  the  isolated  channel  case.  But  a  14cm  spacing  is  well  in  excess 
of  the  “effective  width”  mentioned  earlier,  a  fact  which  would  seem  to  suggest  that  the  array  of  three 
channels,  if  all  operated  together,  would  not  produce  an  appealingly  aperture-efficient  configuration.  But 
that  hypothesis  is  mistaken,  as  shown  in  the  model  in  Figure  15.  The  gain  rises  to  18.72  dB!  Recall  that 
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the  non-interacting  combination  of  three  100%-isolated  channels  of  this  type  should  have  increased  to 
13.244  dB  +  101ogio(3)  =  18.015  dB.  So.  where  did  the  extra  +0.7  dB  come  from?  {continues,  page  14} 


Far  Field:  Gain  Pattam  In  H-plana  (dB) 


Concerns  compared  to  single  isolated  channel: 

•  Performance  seriously  degraded  due  to  coupling. 

•  H-plane  gain  peak  shifted  to  28.0°  (was  30.0°) 

•  Gain  fell  to  9.04  dB  (vs.  13.244  dB  for  the  same 
but  isolated  single-channel  configuration). 

•  Only  ~62%  of  source  power  radiated.  [38%  of  Psrc 
was  variously  reflected  or  lost  into  the  termina¬ 
tions  of  the  driven  and  other  two  channels.] 


Figure  10.  Three  tightly-spaced  channels,  center-channel  driven.  Chan,  spacing:  8  cm. 
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Figure  11.  Selected  curves  from  Figure  4,  showing  falling  S2i  as  channel  separation  rises. 


11 


Gain  Pattern  in  H-plane  (dB) 


|  Ex  |  along  centers  of  each  channel 


Far  Field:  Gain  Pattern  in  H-plane  (dB) 


Compared  to  center-channel  driven  at  8cm  spacing; 

•  Performance  improved  due  to  less  coupling 

*  H-plane  gain  peak  still  at  ~  28.0°  (instead  of  30°) 

*  Gain  is  10.05  dB  (vs.  9.04dB  dB  for  8cm  channel 
spacing)  but  was  13.244dB  for  isolated  chan. 

•  ~78%  of  Psrc  radiated,  (vs  ~62%  for  8cm  channel 
spacing)  but  was  93.8%  for  isolated  chan. 


Figure  12.  Three  channels,  with  only  center-channel  driven.  Chan,  spacing  =  10  cm. 


freq(l)=9.815e8  Slice:  Electric  field,  x  component  (V/m) 


Slice:  Electric  field,  x  component  (V/m)  -0-5 


Ck  tom  3D_FAWSEA_couplinR  C'tion  Ucm  mph 


Far  Field;  Gain  Pattern  In  H-plane  (dB) 


Compared  to  center-channel  driven  at  8cm  spacing: 

•  Performance  improved  due  to  less  coupling 

•  H-plane  gain  peak  is  at  ~28.75°  (instead  of  30°) 

•  Gain  is  11.65  dB  (vs.  9.04dB  dB  for  8cm  channel 
spacing)  but  was  13.244dB  for  isolated  chan. 

•  ~88%  of  Psrc  radiated,  (vs  ~62%  for  8cm  channel 
spacing)  but  was  93.8%  for  isolated  chan. 


Figure  13.  Three  channels,  with  only  center-channel  driven.  Chan,  spacing  =  12  cm. 
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freq(l)=9  815e8  Slice:  Electric  field,  x  component  (V/m) 
Slice:  Electric  field,  x  component  (V/m) 


Slice  view,  Ex 

(in  central  plane) 


P,«./P»<  =  90.68% 

CF.om  5D_ F AWSEA _coopiui t_Cttoa_  1 4cm npk 


Far  Field  Gain  Pattern  in  H-plan*  (dB) 


Figure  14.  Three  channels,  with  only 


Compared  to  center-channel  driven  at  8cm  spacing: 

•  Performance  improved  due  to  less  coupling 

•  H-plane  gain  peak  is  at  ~29.5°  (instead  of  30°) 

•  Gain  is  13.15  dB  (vs.  9.04dB  dB  for  8cm  channel 
spacing).  It  was  13.244dB  for  isolated  chan. 

•  ~90.7%  of  P$rc  radiated,  (vs  ~62%  for  8cm  channel 
spacing).  It  was  93.8%  for  isolated  chan. 


center-channel  driven.  Chan,  spacing  =  14  cm. 
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Far  Field:  Gain  Pattern  in  H-plane  (dB) 


Compared  to  all  channels  driven  at  8cm  spacing: 

•  Performance  improved  due  to  reduced  coupling. 

•  H-plane  gain  peak  now  at  29.5°  (was  at  28.8°  for 
8cm  spacing) 

•  Gain  =  18.72  dB  (vs.  16.48  dB  when  combined 
with  8cm  spacing). 

•  ~86%  of  source  power  radiated. 


Figure  15.  Three  channels,  all  driven.  Channel  spacing  =  14  cm. 
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The  answer  is  that  as  we  increased  the  channel-to-channel  spacing  (and  as  we  observed  the  center  channel 
gain  increasing)  the  gains  associated  with  the  side  channels  were  actually  increasing  faster,  and  actually 
exceeded  that  of  the  isolated-channel  case.  In  other  words,  the  mutual  coupling  of  each  side  channel  to 
the  other  two  channels  enhanced  the  gain.  Though  initially  surprising,  this  behavior  is  hardly  without 
precedent  in  antenna-arrays;  it  is  due  to  constructive  interference  from  parasitic  channels  -  analogous  to 
the  sort  of  behavior  used  to  advantage  in  Yagi-Uda  antennas,  for  example.  Figure  16  reveals  this  effect 
further  by  modeling  the  case  where  the  spacing  is  once  again  14cm,  but  only  one  side-channel  is  driven. 
The  gain  of  this  side  channel  (in  interaction  with  the  two  parasitic  channels  off  to  one  side)  is  14.4dB, 
which  is  about  1.1 6dB  more  than  that  from  a  single  channel  alone. 


freq(l)=9.815e8  Slice:  Electric  field,  x  component  (V/m)  w 

Slice:  Electric  field,  x  component  (V/m)  *0-5 

Slice  view,  Ex 


(driven  side-chan  mid-plane)  \  1 


from  'D_FAWSEA_c«pfc«g_Lcfto._  1  -km  n*>h 


|EX|  along  centers  of  each  channel 


Far  FwM  Gain  Pattern  in  H  piane  (dB) 


Concerns  compared  to  single  isolated  channel: 

•  Coupling  is  evident  to  neighboring  channels. 

•  H-plane  =  14.4  dB  (vs.  13.244  dB  for  an  isolated 
single-channel  configuration).  ->  In  this  case,  the 
coupling  enhanced  the  gain,  for  that  channel. 

•  92.7%  of  source  power  radiated. 


Figure  16.  Three  channels,  only  one  side  channel  driven.  Channel  spacing  =  14  cm. 


Now,  despite  the  arguably-unexpected  18.72  dB  gain  computed  (per  Figure  15)  for  the  all  three  channels- 
driven  design  with  14cm  spacing,  this  design  is  still  somewhat  problematic.  The  total  aperture  width  (if 
we  associate  +4. 5cm  extensions  on  each  side,  per  Weff  noted  earlier)  is  -28+9  =  37cm,  yielding  a  physical 
aperture  area  ~0.74m2.  The  aperture  efficiency  corresponding  to  this  18.72  dB  gain  at  981.5  MHz  for  that 
much  area  follows  immediately,  as  74.7%.  That’s  respectable,  although  we  have  done  somewhat  better  in 
other  FAWSEA  designs.  We  are  more  concerned  about  the  imbalanced  aperture  field  distribution. 


Figure  17  shows  slices  of  IEI  at  2,  4,  and  6  cm  from  the  aperture.  The  magnitude  of  the  field  is  more  non- 
uniform  spatially  than  desirable,  a  condition  which  could  encourage  surface  breakdown. 


▼  0.022  ▼  0.037  ▼  0.053 


+2cm  from  aperture  +4cm  from  aperture  +6  cm  from  aperture 

Figure  17.  |E|  (kV/cm)  in  vicinity  of  the  aperture  for  the  Figure  15  design,  if  Psrc  ~1  GW. 
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Generally,  it  is  desirable  to  ensure  more-uniform  aperture  surface  fields  and  to  reduce  the  imbalance  in 
the  behavior  of  interior  vs  exterior  channels.  For  example,  consider  a  design  like  the  one  in  Figure  15,  but 
with  Rfii  increased  to  4cm  and  the  gap  reduced  to  1cm,  again  yielding  14cm  for  the  channel-to-channel 
separation.  (Strictly  speaking,  we  should  now  re-compute  the  grill  wire  diameters,  but  let’s  just  keep  them 
unchanged  for  now.)  Figure  18  shows  the  predicted  performance  after  this  change,  which  also  happens 
to  deliver  0.5  dB  more  gain  than  the  arrangement  in  Figure  15. 

(No  change  to  grill  wire  diams)  |  Ej  a|ong  centers  of  each  channe| 


0  From  tD  FAWSEA  coupling  alltoo  14cm  tprcuil  mph 


Far  Field:  Gain  Pattern  In  H-plane  (dB) 


With  Rt:i=4.0cm.  gap  =  1cm,  chan  sep  =  14cm: 

•  Performance  improved. 

•  H-plane  gain  peak  at  29.5° 

•  Gain  =  19.22  dB. 

•  ~96%  of  source  power  radiated. 


increased  to  4cm.  Channel  spacing  =  14  cm. 


Figure  18.  Three  channels,  all  driven,  Rfn 


Also,  Figure  19  shows  a  smoother  distribution  of  IEI  near  the  aperture  than  that  appearing  in  Figure  17. 


freq(l)=9.815e8  Multislice 


(kV/cm)  i 

▲  15.44 


+2cm  from  aperture 


+4cm  from  aperture 


+6  cm  from  aperture 


Figure  19.  |E|  (kV/cm)  in  vicinity  of  the  aperture  for  the  Figure  18  design,  if  Psrc  -1  GW. 


We  will  continue  the  discussion  of  mutual  coupling  and  the  challenges  and  opportunities  it  represents  in 
regard  to  the  optimization  of  aperture  field  distributions  in  our  next  report. 
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3.2.  Compensating  for  Modified  Phase-velocities  in  Leaky-wave  Channels 

Our  earlier  analyses  (described  in  our  earlier  reports)  set-aside  the  question  of  explicitly  computing  the 
phase  velocity  of  the  wave  in  a  leaky-guide,  which  differs  somewhat  from  that  in  a  closed-wall  guide. 

For  only  slight  leakage  the  difference  is  small,  but  it  becomes  more  pronounced  if  the  grill  is  leakier.  For 
theoretical  completeness  and  to  speed-up  the  design  process,  the  necessary  adjustment5  is  described  here. 


Recall  that  for  a  uniform  grill  and  wave 
incident  at  angle  0  relative  to  normal, 
we  previously  noted  that  the  complex 
reflection  coefficient  at  the  grill  was 
given  by: 

"2'v2 


R 


■Zo+x: 


Equivalent  PEC  wall,  from 
phase-velocity  perspective 

Leaky  wire  ->  •  • 
grill 


xtra 


•  •  •  •  #T|**  • 

♦  ♦  A  ♦  ♦  ♦  ♦ 

.  ♦  ♦  r  ♦  ♦  ♦  ♦ 

Vi  ♦  ♦  I  ♦  ♦  ♦  ♦ 

♦  1  ♦  ♦  ♦  ♦ 

♦  ♦  ♦  ♦  ♦ 

: 


(z0+j(xa-xj+xy 

where  expressions  for  Xa  and  Xb  are 
given  by  Marcuvitz  (1951)  and  (for  his 
expressions)  Z0  can  be  normalized  to 
unity.  In  fact,  computation  of  Rg  was 
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PEC  back  wall 

Figure  20,  Concept  of  a  wrfi/a/perfectly-conducting 
wall  for  the  phase-velocity  computation.  (Note:  this 
is  an  H-plane  view;  E  is  normal  to  the  page) 


already  included  in  our  previously-documented  custom  MatLab  functions.  The  next  step  to  computing 
the  detailed  impact  upon  the  phase  velocity  within  the  leaky  guide  is  to  recognize  that,  in  this  context ,  the 
grill  behaves  very  much  as  if  it  were  equivalent  to  a  perfectly-conducting  wall  moved  slightly  farther  out 
from  the  opposite  side  wall  of  the  waveguide.  (See  Figure  20  for  the  basic  idea.)  Thus,  in  practice,  all  we 
need  to  do  is  compute  the  value  of  Sxtra .  Once  we  have  that,  we’ll  know  where  to  place  our  grill  wires  to 
match  the  phase  velocity  in  the  feeding  waveguide.  It  is  clear  that  the  distance  Sxtra  corresponds  to  an 
equivalent  phase  factor  =  exp(2j*ko*Sxtra/cos0).  But  since  the  reflection  coefficient  at  a  PEC  wall  is  real , 
the  imaginary  part  of  Rg  (if  its  magnitude  is  normalized)  must  come  solely  from  this  term.  Thus: 


imagy 
I  R„ 


=  sin \2k0Sxtra  / cos  0  ) ,  which  immediately  yields: 


e  cos#  . 
5.na  = - sin 


2  h 


0 


imag[ 

RJ 

Rg 

Of  course,  since  Rg  varies  as  a  function  of  position  if  we  change  the  diameters  (or  spacing)  of  the  wires  as 
a  function  of  position,  then  any  generalized  Sxtra  will  also  vary  as  a  function  of  position,  suggesting  (as  we 
have  previously-noted)  correction  methods  such  as  tilting  the  back-wall  and/or  otherwise  shifting  the  grill 
wire  positions  relative  to  the  back-wall  to  hold  the  phase  velocity  -constant  throughout  the  channel.  For 
future  reference,  a  custom  MatLab  function  to  compute  8xtra  via  the  above  rule  is  provided  in  Figure  21. 


%  function  DeltaXtra  =  DeltaXtra (Angln, a, d, f req) 

%  Angln  =  angle  of  inc .  rel  to  normal,  in  radians 
%  a  =  wire  spacing  in  meters 
%  d  =  wire  diameter  in  meters 
%  freg  =  freguency  in  Hz 

%  Computes  distance  from  grill  to  eguiv  PEC  wall. 

%  Used  for  FAWSEA  design.  E  is  parallel  to  wires. 
%  Prepared  by  R.  Koslover,  SARA,  Inc. 

%  June,  2014 

%  See  Marcuvitz,  Waveguide  Handbook,  p.286. 

%  See  also  Quarterly  Report  #3,  SARA  Proj :  Navy41. 
function  DeltaXtra  =  DeltaXtra (Angln, a, d, freg) 
ct  =  cos (Angln) ; 
st  =  sin (Angln) ; 

c_const  =  2 . 99792458e8 ;  %c  in  m/s 

lambda  =  c_const/f reg;  %  free  space  wavelength 


suml  =  0; 
sum2  =  0; 
for  m  =  1:1:1000 

suml  =  suml+l/sgrt (mA2+2*m*a*st/lambda- 
(a*ct/lambda) A2) -1/m; 

sum2  =  sum2+l/sgrt (mA2-2*m*a*st/lambda- 
(a*ct/lambda) A2) -1/m; 
end 

%  Normalized  to  Z0  =  1  everywhere 
%  Compute  Xa  using  sum  from  Marcuvitz,  p.  286 
Xa  =  (a*ct/lambda) * (log (a/ (pi*d) ) +0 . 5* (suml+sum2) ) ; 
Xb  =  (a*ct/lambda) * (pi*d/a) ^2 ; 

Zeg  =  j* (Xa-Xb) - (Xa^2) / (-1+j* (Xb-Xa) ) ; 
gam  =  (Zeq-1 ) / (Zeq+1 ) ; 

DeltaXtra  = 

(ct*lambda/ (4 . 0*pi) ) *asin (imag (gam) /abs (gam) ) ;  %  meters 
end 


Figure  21.  Listing  of  custom  MatLab  function  ‘DeltaXtra’  that  computes  Sxtra. 


5  This  calculation  is  only  strictly  valid  in  the  ideal  2D  case  (H-plane  analyses  only)  with  no  dielectric  window,  but 
may  be  generalizable,  as  will  be  discussed  shortly. 
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A  plot  of  computed  Sxtra  vs  wire  diameters  at  selected  wire  spacings  is  provided  in  Figure  22. 


Figure  22.  Computed  8xtra  vs  wire  diameter  for  various  wire  spacings, 
scaled  to  the  tree-space  X,  for  0inc=3O°. 

Although  we  have  derived  8xtra  from  the  imaginary  component  of  the  wire-grill  reflection  coefficient  Rg , 
this  could  be  generalized ,  in  principle ,  to  use  the  imaginary  component  of  the  overall  reflection  coeffi¬ 
cient,  the  latter  found  via  the  wave-matrix  formalism  (see  previous  reports),  thus  effectively  accounting 
for  curved  channel  edges  in  the  E-plane,  a  dielectric  window,  and  even  other6  factors  when  positioning 
the  wire  grill.  However,  it  is  unclear  if  additional  adjustments  to  the  grill  or  opposite  wall  that  would 
follow  from  doing  those  more-complete  calculations  would  prove  significant.  We  may  revisit  this  later. 

3.3.  Aperture  Options  for  Forward-traveling  Leaky-wave  HPM  Antennas 

The  far-field  gain  of  any  aperture-type  antenna  has  an  upper-bound  which  is  generally  only  realizable 
with  an  ideal  distribution  of  aperture  illumination.  For  flat  apertures,  the  most  ideal  illumination  is 
100%  uniform6 7  in  terms  of  magnitude,  polarization,  and  phase.  Traveling- wave  leaky- wave  antennas 
cannot  deliver  uniform  phase,  but  they  can  deliver  fairly-uniform  magnitude  and  polarization,  with 
nearly-linear  phase  along  the  aperture.  The  latter  results  in  a  non-zero  tilt  angle  of  the  beam  relative  to 
the  aperture  normal.  Considered  separately  from  other  factors,  this  imposes  a  factor  of  cos(4>tiit)  on  the 
numerical  gain,  where  <\>tnt  is  the  beam  tilt  angle.  This  reduces  the  effective  area:  Aleaky  =  Ageo  cos  (j)tih . 

Source  bandwidths,  cutoff  considerations,  and  finite-duration  pulse  considerations  typically  lead  us  to 
choose  (|)tiit  >  -20°  or  so  for  practical  HPM-capable  designs.  In  particular,  selecting  c|)tilt  -25-35°  often 
provides  a  satisfying  tradeoff  between  bandwidth  and  gain.  Note  that  within  that  range,  the  tilt-related 
efficiency-factor  ranges  from  cos(25°)  =  90.6%  to  cos(35°)  =  81.9%.  However,  there  is  more  to 
describing  antenna  patterns  than  noting  the  peak  gains  achievable,  and  there  are  many  more  apertures  of 
interest  than  those  that  are  simple  and  flat. 


6  For  example,  one  could  potentially  account  for  the  curved  section  between  the  grill  and  aperture  in  a  RAWSEA. 

7  We  exclude  “super-gain”  class  antennas,  as  their  extremely  high-Q  configurations  preclude  HPM/DEW  appli¬ 
cations  and,  in  fact,  most  other  applications  as  well.  Note  also  that  the  far-fields  of  non-aperture  (e.g.,  wire-type) 
antennas  can  also  be  characterized  by  the  approach  to  be  described  here,  but  it  is  generally  less  useful  to  do  so  than 
to  employ  alternative  methods  integrating  over  their  surface  currents. 
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We  are  currently  cataloging  some  of  the  aperture  shapes  we  consider  to  be  of  most  interest,  along  with 
associated  antenna  patterns  and  their  anticipated  practically-achievable  aperture  efficiencies.  Of  course, 
curved  apertures  are  usually  less  aperture  efficient  than  flat  ones,  but  may  still  deliver  quite-respectable 
and  operationally-practical  gain,  especially  when  compensating  feed  schemes  are  introduced  (such  as  we 
have  noted  in  earlier  analyses).  For  reference,  several  of  the  aperture  geometries  in  our  current  study  are 
shown  in  Figure  23,  along  with  example  representative  aperture  fields,  although  the  latter  (in  these 
particular  examples)  are  shown  without  distributed  phase  delays  or  other  compensation  for  gain-sapping 
aperture  curvature.  Note  also  that  the  sizes  and  aspect  ratios  of  the  apertures  shown  here  were  chosen 
arbitrarily,  i.e.,  they  are  not  intended  to  represent  any  particular  candidate  HPM  DEW  platforms. 


H-plane  curvature  with  a  pinch 


Both  E  &  H  curvature  with  a  pinch 


Figure  23.  Some  of  the  leaky-wave  aperture  geometries  currently  under  investigation. 
(Note:  Only  the  apertures  are  highlighted  above,  for  simplicity.) 
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We  are  exploring  the  far-field  patterns  that  can  be  generated  by  forward-traveling,  leaky-wave  field 
distributions  in  apertures  that  are  variously  flat,  curved  in  the  E-plane,  curved  in  the  H-plane,  curved  in 
both  planes,  or  curved  and/or  pinched  more  creatively.  Field  re-distribution  schemes  to  compensate,  at 
least  in  part,  for  the  gain-degrading  effects  of  aperture  curvatures  are  also  being  studied.  For  future 
reference,  we  note  that  for  computation  of  the  far-fields  from  all  these  antennas,  we  are  currently  using 
the  following  balanced  aperture-integration  approximation  formula8,  which  is  among  the  more 
practical/useful9  expressions  available: 

EP  =  ^ K  x  jj  [n  x  Ea  -  r/rp  x  (n  x  Ha  )]exp  (  jkra  ■  rp  )dS 

s 

where  Ep  is  the  computed  electric  field  at  the  far-field  point  p ,  Ea  and  H a  are  the  electric  and  magnetic 

complex  vector  fields  on  the  aperture  (or  surface  S  enclosing  the  antenna),  h  is  the  unit  outward  normal 
on  the  aperture  (or  surface  5),  77  is  the  impedance  of  free  space  (-377  Ohms),  k  is  the  free  space  wave 
number  (2nlX),  ra  is  the  vector  from  the  origin  to  the  aperture  field  point,  rp  is  the  vector  from  the  origin 

to  the  radiated  field  point,  with  rp  =  rp  /\rp  |  =  a  unit  vector  along  rp  .  In  cartesian  coordinates, 

rp  =  x  sin  6p  cos  (/>p  +  y  sin  6p  sin  </>  +  z  cos  6p  .  [Note:  for  a  uniform-field  flat  aperture  in  the  xy-plane, 
the  above  surface  integral  simplifies  greatly  and  it  is  fairly  easy  to  show  that  the  resulting  expression  for 
Ev  at  9P  =0  is  consistent  with  on-axis  gain  Gideal  =  A7lAgeol  .]  Of  course,  analytic  or  semi-analytic 

expressions  for  E-plane  and  H-plane  cuts  of  far-field  patterns  for  a  variety  of  idealized  aperture  antenna 
distributions  can  be  found  in  the  literature10,  although  most  are  limited  to  planar  apertures,  especially 
those  with  rectangular  or  circular  boundaries.  (Recall  also  the  previously-noted  the  works  of  Honey, 
Ishimaru,  and  others  in  analyzing  leaky-wave  antennas  and  apertures  conformal  to  ideal  flat  or  cylindrical 
surfaces.)  For  greater  generality,  our  approach  here  is,  and  will  continue  to  be,  mostly  numerical.  We 
have  previously  implemented  the  aforementioned  aperture-field  integration  expression  for  Ev  in  a  variety 
of  circumstances  and  in  some  cases  have  used  more  complicated  forms  to  compute  radiating  near-fields 
and  time-domain  responses,  although  those  are  not  expected  to  be  needed  here. 

In  our  next  report,  we  plan  to  include  a  summary  of  the  flat  and  curved  aperture  shapes  studied,  more 
information  about  the  leaky-waveguide  structures  that  are  so  essential  to  fitting  these  apertures  and 
enabling  their  practical  realization,  and  their  predicted/anticipated  performance  characteristics  of  these 
novel  antennas,  including  (but  not  limited  to)  their  predicted  far-field  patterns. 


4.  DISCUSSION,  CONCLUSIONS,  AND  RECOMMENDATIONS 

We  are  pleased  to  report  that  the  work  performed  during  this  3rd  quarter  of  the  R&D  program  has 
continued  to  advance  this  important  technology  and  to  extend  the  tool  set  available  to  would-be  designers 
of  leaky-wave  based  HPM  antennas.  We  look  forward  to  both  extending  the  theory  and  developing  and 
documenting  in  our  next  report  a  set  of  representative  example  designs  to  better  guide  the  engineering  of 
these  antennas.  We  appreciate  ONR’s  continued  support  for  this  R&D. 


8  S.  Silver,  Microwave  Antenna  Theory  and  Design.  See  especially  Sec.  5-11,  “The  Aperture  Field  Method,”  and 
Sec.  5-12,  “The  Fraunhofer  Region.”  From  1st  Ed,  publ.  by  office  of  Scientific  Research  and  Development, 
National  Defense  Research  Committee,  NY,  1949.  See  http ://www. \ lab.org/ir/MITSeries/V  1 2 .PDF  . 

9  This  particular  expression  is  especially  convenient  numerically,  since  it  contains  no  derivatives. 

10  See,  among  others,  Jull,  E.V.,  “Radiation  from  Apertures,”  Chap.  5  of  Antenna  Handbook:  Theory,  Applications, 
and  Design,  Ed.  by  Y.T.  Lo  and  S.W.  Lee,  Van  Nostrand  Reinhold,  NY,  1988. 
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